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The Arnowitt-Deser-Misner (ADM) evolution equations for the induced metric and the 
extrinsic-curvature tensor of the spacelike surfaces which foliate the space-time manifold 
in canonical general relativity are a first-order system of quasi-linear partial differential 
equations, supplemented by the constraint equations. Such equations are here mapped 
into another first-order system. In particular, an evolution equation for the trace of the 
extrinsic-curvature tensor K is obtained whose solution is related to a discrete spectral 
resolution of a three-dimensional elliptic operator V of Laplace type. Interestingly, all 
nonlinearities of the original equations give rise to the potential term in V. An example of 
this construction is given in the case of a closed Friedmann-Lemaitre-Robertson- Walker 
universe. Eventually, the ADM equations are re-expressed as a coupled first-order system 
for the induced metric and the trace-free part of K. Such a system is written in a form 
which clarifies how a set of first-order differential operators and their inverses, jointly with 
spectral resolutions of operators of Laplace type, contribute to solving, at least in principle, 
the original ADM system. 
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The canonical formulation of general relativity relies on the assumption that space- 
time (M, g) is topologically E x R and can be foliated by a family of spacelike surfaces E t , 
all diffeomorphic to the three-manifold E. The space-time metric g is then locally cast in 
the form 

g = — (jV 2 - NiN*)dt ®dt + N i (dx i <g> dt + dt <g> da;*) + h ij dx i <g> ofe j , (1) 

where iV is the lapse function and N l are components of the shift vector of the foliation 
[1-3]. The induced metric hij on E t and the associated extrinsic-curvature tensor Kij turn 
out to obey the first-order equations [1-3] 

= -2NK it + N ib + N lh , (2) 



^R ij + K ij (tvK)-2K im K, 



in 



+ 



N m K tjlm + N^K jm + N^K im 



(3) 



where the stroke denotes covariant differentiation with respect to the induced connection 
on E*. Moreover, the constraint equations hold. For Einstein theory in vacuum they read 

K^-itvK^^O, (4) 

(3) i? - KijK ij + (trK) 2 w 0, (5) 

where ~ is the weak-equality symbol introduced by Dirac to denote equations which only 
hold on the constraint surface [3,4]. 

Equation (3) is quasi-linear in that it contains terms quadratic in the extrinsic- 
curvature tensor, i.e. 

N(K ij (trK)-2K im K J m ). 

We are now aiming to obtain from eqs. (2) and (3) another set of quasi-linear first-order 
equations. For this purpose, we contract Eq. (2) with and Eq. (3) with h % K This 
leads to (hereafter —A = — ^ '* is the Laplacian on E t ) 

O 7 

RlJ ~dt L = - 2NK v Rij + 2RijN i\j, ( 6 ) 
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r)K 

fcvZfLy =-AN + N 



dt 



(3) 



R+(trK) 2 -2K im K i 



+ 



N m {trK)\ m + 2N^K i m 



(7) 



It is now possible to obtain an evolution equation for (trK) , because 



dV 1 dt 13 dt 



(8) 



The second term on the right-hand side of Eq. (8) is given by Eq. (7), whereas the first 
term is obtained after using the identity 



d_ 
di 



d 

di ( 



0, 



which implies 



and hence 



dklJ =-h^ dhpl 



dt 

dh ij 



dt ' 



(9) 



(10) 



dh- 

dt KlJ ~ dt K ■ 



(11) 



By virtue of Eqs. (6)-(8) and (11), and imposing the constraint equation (5), we get 



d 



-(trK) = - (A - KyK**) N + iV m (trK), m , 



(12a) 



which can also be cast in the form (here V m = i m ) 



d 



- - N m V m ) (trK) = (- A +K ij K i *) N. 



(126) 



This is a first non-trivial result because it tells us that, given the operator of Laplace type 
[5] 

V = - A+K i:j K ij , (13) 

the right-hand side of Eq. (12b) is determined by a discrete spectral resolution of V. By 
this one means a complete orthonormal set of eigenfunctions belonging to the eigenvalue 
A, so that, for each fixed value of t, the lapse can be expanded in the form 



N(x,t) = J2Cxf p x (x,t), 



(14) 



with Fourier coefficients C\ given by the scalar product 



C A =(/£,iv). (15a) 

This point is simple but non-trivial: for each fixed value of £, the restriction of the lapse 
to E f becomes a function on T, t only, and hence the Fourier coefficients in (15a) read 

Cx = [ (f p x TNVhd 3 x = C x ,t, (156) 

where the star denotes complex conjugation. In the application to the initial- value problem, 
one studies the lapse at different values of t, and hence it is better to write its expansion 
in the form (14), with the time parameter explicitly included. The integration measure for 
C\, however, remains the invariant integration measure on S f . 

The existence of discrete spectral resolutions of V is guaranteed if E t is a compact 
Riemannian manifold without boundary [5], which is what we assume hereafter. Thus, Eq. 
(12b) can be re-expressed in the form 

ijt ~ Nm ^ m ) tri ^) = E AC v /f(o?,t). (16) 



On denoting by L the operator 



L = -^-N m V m , (17) 



and writing G(x, t; y, t') for the Green function of L, Eq. (16) can be therefore solved for 
(trK) in the form 

tvK(x,t)= [ G(x,t;y 7 t')J2*C x , t ,f p x (y,t')dv(y,t'), (18) 

JT t i xR A 

where d/j,(y, t') is the invariant integration measure on T, t / xR. Note that the eigenfunctions 
f x are different at different times and hence, to compute them at all times, one has to 
integrate the equations of motion so that the metric and the extrinsic-curvature tensor 
are known at all times. Thus, Eq. (18) does not reduce the amount of entanglement of 
the original equations (2) and (3). A relevant particular case is obtained on considering 
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the canonical form of the foliation [6] , for which the shift vector vanishes. The operator L 

d_ 



reduces then to 4t, and its Green function can be chosen to be of the form 



G c {t,t')= 1 -[e{t-t')-e{t'-t)\ (19) 

where 6 is the step function. Such a Green function equals \ for t > t' and — \ for 
t < t', and hence its "jump" at t' equals 1, as it should be from the general theory. It is 
sometimes considered in a completely different branch of physics, i.e. the theory of solitons 
and nonlinear evolution equations. 

A cosmological example shows how the right-hand side of Eq. (16) can be worked 
out explicitly in some cases. We are here concerned with closed Friedmann-Lemaitre- 
Robertson-Walker models, for which the three-manifold E reduces to a three-sphere of 
radius a. It is indeed well known that the space of functions on the three-sphere can be 
decomposed by using the invariant subspaces corresponding to irreducible representations 
of 0(4), the orthogonal group in four dimensions. These invariant subspaces are spanned 
by the hyperspherical harmonics [7], i.e. generalizations to the three-sphere of the famil- 
iar spherical harmonics. The scalar hyperspherical harmonics Q( n \xi ®-> V 9 ) form a basis 
spanning the invariant subspace labeled by the integer n > 1 corresponding to the n-th 
scalar representation of 0(4). The label (n) of refers both to the order n and to other 
labels denoting the different elements spanning the subspace. The number of elements 
is determined by the dimension of the corresponding 0(4) representation. On using 
the general formula for 0(2k) representations, the dimension of the irreducible scalar rep- 
resentations is found to be d,Q(n) = n 2 . The Laplacian on a unit three-sphere when acting 
on 0/ n ) gives 

-AQ (n ) = (n 2 -l)Q (n) n=l,2,.... (20) 

Any arbitrary function F on the three-sphere can be expanded in terms of these hyper- 
spherical harmonics as they form a complete set: 

F = J2^)Q {n) > (21) 

(n) 

where the sum (n) runs over n = 1 to oo and over the n 2 elements in each invariant 
subspace, and the are constant coefficients. The lapse function on the three-sphere is 
therefore expanded according to (cf. (14)) 

N = J2 c (n)Q (n \ (22) 

(n) 
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and bearing in mind that KijK^ = 3 on a unit three-sphere we find, on a three-sphere of 
radius a, 




(n) 



where (here a = a(t)) 




n = 1,2, ... . 



(24) 



We have therefore found that all nonlinearities of the original equations give rise to the 
potential term in the operator V defined in Eq. (13), which is a positive-definite operator 
of Laplace type (unlike the Laplacian, which is bounded from below but has also a zero 
eigenvalue when n = 1). 

The interest in the evolution equation for the trace of the extrinsic-curvature tensor 
is motivated by a careful analysis of the variational problem in general relativity. More 
precisely, the "cosmological action" is sometimes considered, which is the form of the 
action functional I appropriate for the case when (trK) and the conformal three-metric 
are fixed on the boundary. One then finds in c = 1 units [8] 



and this formula finds important applications also to the quantum cosmology of a closed 
universe, giving rise to the if-representation of the wave function of the universe [9]. 

The program we have outlined in our letter shows an intriguing link between the 
Cauchy problem in general relativity on the one hand [10-13], and the spectral theory of 
elliptic operators of Laplace type on the other hand [5], which is obtained by exploiting 
the nonlinear form of the ADM equations and the nonpolynomial form of the Hamiltonian 
constraint (5). On denoting by <7y the trace- free part of the extrinsic-curvature tensor: 



we therefore obtain a set of equations, equivalent to (2) and (3), in the form (here L 
denotes the inverse of the operator L defined in (17)) 




(25) 



(Tij = Kij - -hij(trK) 



(26) 



trK = L~ x VN, 



(27) 



dha 



-Nhij(tTK) - 2Naij + 2N m , 



(28) 



dt 



6 



da 



>3 



1 



0f n tJ N - -N m (tvK) + N m V m cr t j + 2JV$ \a j)m + -h j)m (tvK) 
where the operator of Laplace type fiy is given by 

Q.. = -Vj-Vi - hi l3 V + ^(trK) 
+ ^h tJ (trK) 2 + WRij - 2a irn a™. 



(29) 



(30) 



It now appears desirable to understand to which extent the differential-operator view- 
point plays a role in solving the coupled system (27)-(30). For this purpose, we consider a 
relevant particular case, i.e. the canonical foliation studied in Ref. [6], for which the shift 
vector vanishes. The operator L reduces then to and, on defining the operators 



Q = ^- + = l+ \nl~ 1 vn, 

at 3 3 



(31) 



one finds the equations 



dt 3 



L-^L-'VN, 



A 



■ij 



(32) 
(33) 



S H — (L~ 1 VN) 2 Q~ 1 N 



(34) 



h 



-2Q" 1 iV(j i? -, 



(35) 



supplemented, of course, by Eq. (27). In other words, once Eq. (34) is solved for the trace- 
free part of the extrinsic-curvature tensor of Et, the induced metric on is obtained from 
Eq. (35). This form of the coupled first-order ADM system of equations has possibly the 
merit of stressing the need to invert the operators L, Q and S to find a solution for given 
initial conditions. Approximate solutions, to the desired accuracy, will correspond to the 
construction of approximate inverse operators Q~ x and S~ x . For this purpose, it can 
be useful to re-express Eq. (34) in the form 



I + -S^iL^VNfQ^N 
3 



(34') 
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The action of the inverse operator L _1 on VN is already given by Eq. (18), but the inverses 
Q~ x and S~ x make it necessary to develop an algorithm for their exact or approximate 
evaluation. More precisely, one has 

Q = LU+ 2 -L-'NL-^N) , (31') 



(32') 



and hence 



Q' 1 = (j+ ^L^NL^VN^j L~\ (36) 

S' 1 = (j- ^L^NL^VN^j L~\ (37) 

by virtue of the operator identity (AS) -1 = B~ 1 A~ 1 . From this point of view, two levels 
of approximation seem to emerge: 

(i) The degree of accuracy in the evaluation of Q~ l and S~ x . It is clear from Eqs. (36) 
and (37) that this reduces to finding 

(I + pL^NL^VN)' 1 

where p = | or — g. A series expansion will be, in general, only of formal value. However, 
on using the norms defined in section 7.4 of Ref. [14] , which relies on well known properties 
of Sobolev spaces, if the operator 

T p = pL^NL^VN (38) 

is a bounded operator on a Banach space X with norm ||T p || < 1, the inverse of I + Tp 
exists and the series ^^ =0 (— l) n (T p ) n converges uniformly to (/ + T p ) _1 with respect to 
the norm on the set of bounded maps from X into X, and one can write 

a+^r^E^rra" ( 39 ) 

n=0 

Equation (39) can be therefore used if L~ l is such that T p is a bounded operator with 
\\T P \\ < 1. 
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(ii) The way in which the non-linear term a"£ m cr- m is dealt with. For example, one may 
regard the right-hand side of Eq. (34') as the "known term", and hence consider the 
equation 

a l3 dx l <g> dx 3 - ^F^S' 1 [(Q~ 1 Na ij )VN] dx 1 <g> dx J 

- F- 1 S- 1 {A, 3 N)dx i <g> dx 3 = -2F" 1 ^" 1 (Naima™) dx 1 <g> dx 3 \ (40) 
where F _1 is the inverse of the operator 

F = I + -S~ 1 (L~ 1 VN) 2 Q~ 1 N. (41) 
3 

Equation (40) is an integral equation for cr^-, for which a perturbation approach might be 
useful if one takes 

-2F- 1 ( S" 1 {Na im a 3 m ) dx 1 <g> dx 3 

as the known term mentioned before. Note that we have resorted to the use of the tensor 
product dx 1 <g> dx 3 because in the resulting equation (40) one has a well defined integration 
of a "symmetric two-form" over the space-time manifold, here taken to be diffeomorphic 
to E f x R. Strictly, also Eq. (35) should be written as 

hijdx 1 <g> dx 3 = -2(Q- 1 No lJ ^jdx l ® dx 3 , (35') 

and in all such equations only the symmetric part of the tensor product survives, because 
both hij and o~ij are symmetric rank-two tensor fields. 

Note once more that it would be wrong to regard Eq. (27) as independent of the 
solution of Eqs. (28) and (29), because L~ x VN is only known at all times when the time 
evolution of and Kij has been determined for given initial conditions. Our equations 
(35')-(37) and (40), although not obviously more powerful than previous schemes, prepare 
the ground for an operator approach to the ADM equations for general relativity, and 
hence might contribute to the understanding of structural properties of general relativity. 
In particular, our way of writing the coupled system of non-linear evolution equations might 
lead to a better understanding of the interplay between elliptic operators on the spacelike 
surfaces E t and hyperbolic equations on the space-time manifold. As far as we can see, 
this expectation is supported by the closed Friedmann-Lemaitre-Robertson- Walker model 
discussed before (where cr^- vanishes) , and further examples of cosmological interest might 
be found, e.g. Bianchi IX models describing a closed but anisotropic universe. For open 
universes or asymptotically flat space-times, however, we are not aware of theorems that 
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make it possible to expand the lapse as in (14). In such cases, the counterpart of the 
elliptic theory on Ef advocated in our paper is therefore another open problem. 
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